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Abstract. In this paper we study the problem of extending means to means of 
higher order. We show how higher order means can be inductively defined and 
established in general metric spaces, in particular, in convex metric spaces. As a 
particular application, we consider the positive operators on a Hilbert space under 
the Thompson metric and show that the operator logarithmic mean admits exten- 
sions of all higher orders, thus providing a positive solution to a problem of Petz 
and Temesi [17]. 

1. Introduction 

> 

A mean of order n or n-mean for short on a set X is a function u : X n — > X 
G\ ■ ^ 

satisfying fj,(x, for all x G X. Growing out of foundational papers on 

(N \ 

the subject such as that of Kubo and Ando [TT], the theory of matrix and operator 
\Q • 

means of order 2 has seen a substantial development. However, no such developed 
J ■ theory has arisen for means of higher order. Indeed there is no obvious definition for 
the extension of a mean to the next higher order, although for various specific means 
one can find a variety of candidates that have been put forward. 

Our purpose in this paper is to study a method of extending means to higher 
orders that appears to offer a viable general approach. For the definition and method 
of extending, we favor a general version of the recent approach of Horwitz [9] for 
the case of means on the positive reals and of Ando, Li, and Mathias [1] for the 
case of the geometric mean on the positive (semi) definite Hermitian matrices. This 
approach has also been adopted and generalized beyond the geometric mean by Petz 
and Temesi in [17], |16] , although they only obtain existence for ordered tuples. In 
this paper we show that this approach can be generalized to means on metric spaces 
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and develop the theory of these extensions in this context. The main theorems give 
rather general conditions that guarantee that extensions of all order exist. Our main 
applications involve higher order means of positive operators on a Hilbert space. We 
are particularly interested in those cases in which one starts with a mean of two 
variables and inductively extends it to all dimensions greater than two. 

In Section 2 we present our approach to extending means via limits of the "barycen- 
tric operator." (This method is called "symmetrization" in p2].) Section 3 contains 
a major result of the paper: means that are nonexpansive and coordinatewise con- 
tractive admit extensions to all higher orders. In Section 4 a special case of such 
means are considered, namely convex means for which the mean is a metrically con- 
vex function assigning to any two points a metric midpoint. As explained in Section 
5, Hadamard spaces (for which the metric satisfies the semiparallelogram law) form 
an important class of examples of metric spaces with associated convex mean. 

In Section 6 we develop machinery for showing that certain types of iterated means 
are nonexpansive and coordinatewise contractive. Since many important means arise 
in this fashion, this is a useful and important result. Section 7 presents categorical as- 
pects of means and their extensions. In Section 8 a reverse construction is considered: 
given a mean, when is it an extension of a lower order mean? Connections between 
means, their extensions, and order are developed in Section 9. The paper closes in 
Sections 10 and 11 with the study of means on the space of positive operators on a 
Hilbert space. It is shown that certain iterated means are nonexpansive and coordi- 
natewise contractive, hence extend to higher orders. This is true, for example, of the 
logarithmic mean, as we show in Section 11, and hence it has higher order extensions, 
a conclusion that provides a positive solution of a problem of Petz and Temesi in |17j . 

2. Mean Extensions 

An n-mean on a set X consists of an n-ary operation (function) /i : X n — > X that 
satisfies a generalized idempotency law: fi(x,x, . . . ,x) = x for all x G X. A mean 
is just an n-mean for some n > 2. The mean is symmetric if it is invariant under 
permutations: 

//(x^i), . . . ,x n (n)) = n(xi, ■ ■ ■ ,x n ) for any permutation 7T on {1, . . . , n}. 
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(Note that for 2-means this means that the binary operation given by the mean 
is commutative.) A topological n-mean consists of a Hausdorff topological space 
equipped with a continuous n-mean u. We also call the operation u a topological 
n-mean. 

A principal goal of our work is to "extend" a (symmetric) n-mean on X to a 
(symmetric) n + 1-mean. As mentioned in the introduction we pattern our approach 
after [9] and flj, an approach that has some relation to the process of compounding 
three means in three variables to obtain another such mean (see pf]). 

Definition 2.1. Given a set X and a fc-mean fi : X k — > X, the barycentric operator 
(3 = ^: X k+1 X k+1 is defined by 

/?(x) := (/i(vr^!x), . . . , /i(7r^ fc+1 x)) , 

where x = (x\, . . . , Xk+i) G X k+l and 7T^yX := (x\, . . . , Xj_ 1: Xj+i, . . . , Xk+i) G X k . 

For a topological fc-mean, we say that the barycentric map j3 is power convergent 
if for each x G X k+1 , we have lim n /3 n (x) = (x*, . . . , x*) for some x* G X. 

As a motivating geometric example for the terminology consider the 3-mean in IR 3 
that assigns to any three points the centroid of the triangle for which they are the 
vertices, i.e., the point where the three medians meet. If we take now the four vertices 
of a 3-simplex or tetrahedron in IR 3 , the barycentric operator applied to the 4-tuple 
consisting of the four vertices replaces each vertex with the centroid (barycenter) of 
the face opposite it, the face with vertices the remaining three vertices. Thus one 
may envision the tetrahedron with vertices the four centroids of the four faces as the 
result. Repeating this process, one obtains a shrinking family of tetrahedra whose 
intersection is the barycenter of the original tetrahedron, represented by the 4-tuple 
with all entries equal to that point. 

Remark 2.2. There is an alternative way that the barycentric map of a fc-mean 
may be defined, namely instead of defining the z -coordinate of /3(x) by deleting the 
^-coordinate of x G X k+1 and applying fi, we delete the coordinate i* := k + 2 — i 
and then apply \i. This means that we begin (from left to right) by first deleting 
coordinate k + 1, k, down to 1, instead of beginning by deleting coordinate 1 and 
continuing up to k + 1. We denote this alternative barycentric map by (3*. One 
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may define /3*(x) alternatively by reversing the {k + l)-tuple /3(x). Note, as long as 
the mean /i is symmetric, both methods power converge to the same limit, provided 
one of them power converges. This equality of limit does not hold in general for 
nonsymmetric means however. The theories for (3 and (3* run parallel, so we restrict 
our attention to (3 with a few brief remarks concerning (3*. 

We define our first notion of an extension in terms of the barycentric operator. 

Definition 2.3. A mean v : X k+1 — > X is a (3-invariant extension of fi : X k — > X if 
v o f3fj, = v, that is, 

(2.1) 2/(x) = I/(ji(7T^iX), . . . , /i(7T^ fc+1 x)) for all X = (xx, . . . , x fc+1 ) G X k+1 . 

The notion of a /3-invariant extension was introduced by Horwitz [9] , who called it 
type I invariance. 

Proposition 2.4. Assume that fx : X k —>■ X is a topological k-mean and that the 
corresponding barycentric operator [3 is power convergent. Define fx : X k+1 — >• X by 
/2(x) = x* where lim n /3 n (x) = (x*, . . . ,x*). 

(i) // : X k+1 ^ X is a (k + l)-mean on X that is a (3-invariant extension of \i. 

(ii) Any continuous mean on X k+1 that is a (3-invariant extension of fi must equal 

A- 

(iii) If n is symmetric, so is jl. 

Proof, (i) For x e X, x = (x, . . . , x) £ X k+1 , we have /3(x) = x by the idempotency 
of /i. Thus x = lim n /5 n (x) and hence /t(x) = x, i.e., /i is a mean. Further, we have 

/i(/?(x)) = 7n(lim/3 n (/3(x))) = 7n(lim/3 n+1 (x)) = /i(x), 

n n 

where n 1 is projection into the first coordinate. Thus fx o (3 = ft, i.e., fx defines a 
/3-invariant extension of \x. 

(ii) Suppose that v is a continuous + l)-mean of X that is a /3-invariant extension 
of /i. Since v = vo(3 = vo(3 n (by repeated application of the first equality), for 
x G X k+ \ 

Kx) = !/(^(x)) = z/(/3"(x)) = z/(x*, . . . , x*) = x* = £(x), 
where (x*, . . . , x*) = lim^oo /3 n (x). 
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(iii) If fi is symmetric, then f3 commutes with any permutation applied to the entries 
of x 6 X k+l , hence also /3 n , and thus one obtains the same limit with constant entry 
jl(x) in either case. Hence ft : X h+1 — > X is symmetric. □ 

Remark 2.5. If /3* is power convergent, then a /3*-invariant extenstion is defined in 
the manner of the previous proposition and the analogous proposition holds for (3* . 
However, /3-invariant extensions need not be /3*-invariant and vice-versa. But both 
notions collapse to the same one for symmetric means. 

We seek a notion of mean extension that allows one both to deduce readily that a 
large number of properties transfer from a mean to its extension and also is applicable 
to a wide variety of means. The preceding proposition provides the ingredients for 
this definition. 

Definition 2.6. A (k + l)-mean v is a (3- extension of a topological /c-mean li (or 
^-extends li) if for each x e X k+1 , lim n /3 n (x) = (V(x), . . . , f(x)). In this case we say 
that f3 power converges to v, written /3™ — > v. 

We restate parts of Proposition 12.41 in terms of this definition. 

Corollary 2.7. If power converges, where ii is a topological mean, then it converges 
to a (k + l)-mean ft, which (by definition) is a ^-extension of \i. Furthermore, fi, if 
continuous, is the unique ^-invariant extension of \i. 

Remark 2.8. A. Horwitz [9] and later D. Petz and R. Temesi [17J consider means 
on the positive reals and show that any continuous symmetric 2-mean that is strict 
(min(a, b) < fi(a, b) < max(a, b) for a ^ b) and order-preserving in each variable has a 
power convergent barycentric map, and hence has a unique /3-extension to a 3-mean. 
Petz and Temesi point out that the argument for power convergence extends to higher 
order variables, and thus one can inductively define /3-extensions for all n > 2 [TTJ 
Section 5]. For the arithmetic, geometric, and harmonic means the extensions yield 
the usual corresponding means of n- variables. To check this, one has only to note that 
they are continuous and are /^-invariant extensions, then apply the previous corollary. 
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3. Power Convergence 



In this section we consider properties preserved by /3-extensions and develop suffi- 
cient conditions for a topological mean to (recursively) admit a /^-extension. 

Given X equipped with a /c-mean /i, a subset C is convex if /i(xi, . . . ,Xk) G C 
whenever x±, . . . , G C. 

Lemma 3.1. If a topological mean admits a ^-extension, then any closed set that is 
convex with respect to the mean is convex with respect to the extension. 

Proof. Let fj, : X k — > X be the given mean, and let x±, . . . , Xfc+i G A, a closed ji- 
convex set. Set x := (x±, . . . ,Xk+i). Then by convexity each coordinate of /3(x) is 
in A and by induction each coordinate of /3™(x) is in A. Since A is closed, it follows 
that the coordinate limits, which are all p,(xi, . . . ,Xk+i), belong to A, where ft, is the 
/^-extension. □ 

Recall that the convex hull of a set A is the smallest convex subset containing 
A, and can be obtained by intersecting all convex sets containing A. In a similar 
fashion in the case of a topological mean \x the closed convex hull can be obtained by 
intersecting all closed convex sets containing A or, as follows from the continuity of 
/i, by closing up the convex hull. 

Definition 3.2. A topological mean is locally convex if there exists at each point a 
basis of (not necessarily open) neighborhoods that are convex. A metric topological 
mean is uniformly locally convex if for each e > 0, there exists 5 > such that 
the diameter of the convex hull of A is less than e whenever the diameter of A 
is less than 5. A metric topological mean is closed ball convex if all closed balls 
B £ (x) := {y G X : d(x, y) < e} are convex for all x G X. 

Remark 3.3. Since in a topological mean the closure of a set A is convex whenever 
A is and has the same diameter, the convex sets in the definition of locally convex 
and uniformly locally convex may be taken to be closed if the space is regular, which 
is the case in the metric setting. 

Lemma 3.4. Given a topological mean on a metric space, closed ball convexity implies 
uniform local convexity, which in turn implies local convexity. 
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Proof. Let e > 0. Choose 6 = e/4. Then any set A of diameter less than 5 is 
contained in a closed ball of radius less than 26 around any point of A, which in turn 
has diameter less than 45 = e. If X is closed ball convex, then this closed ball is 
convex and hence contains the closed convex hull of A. The proof that uniformly 
locally convex implies locally convex is straightforward. □ 

The next lemma is an immediate consequence of Lemma 13.11 

Lemma 3.5. Let X be a topological mean with a (3 -extension. If X is locally convex 
resp. metric and uniformly locally convex resp. metric and closed ball convex with 
respect to the given mean, then it is with respect to the extension. 

Lemma 3.6. If v is a (3-extension of the topological mean /x and if X is locally convex 
and regular, then v is continuous. 

Proof. Let /x : X k — > X and let v : X k+1 — > X be the /3-extension. Let x = 
(xx, . . . , Xk+i) G X k+1 , let x* = ^(x), and let U be an open set containing x*. Pick 
a closed convex neighborhood V of x* such that V C U. Since by hypothesis the 
sequence /3 n (x) converges to the k + 1-string with entries x*, we have /3 n (x) G V k+1 
for some n large enough. By continuity of fi and hence of (3 n , there exists W open in 
X k+1 containing x such that j3 n (W) C V k+1 . For any y G W, we have /3 n (y) G V fc+1 , 
and hence /3 m (y) G V fc+1 for all m > n since is convex. Since V is closed it follows 
that u(y) G V. Thus z/ is continuous. □ 

Definition 3.7. Let \x : X k — > X be fc-mean on a metric space X. For x = 
(x±, . . . ,Xk+i) G X fc+1 , set |x| = {x\, . . . ,Xk+i}, the underlying set of the (k + 1)- 
tuple, and define the diameter A(x) of x by 

A(x) = diam|x| = sup{d(xi, xA : 1 < i, j < k + 1}. 

The mean /x is weakly (3 -contractive if for each x G X fc+1 , we have lim n A(/9 n (x)) = 0. 
For < p < 1, we say that /x is coordinatewise p- contractive if for any x, y G X k that 
differ only in one coordinate, say Xj ^ yj, 

d(/x(x),/x(y)) < pd(xj,yj) 

Lemma 3.8. If \x : X fc X is a coordinatewise p-contractive mean for < p < 1, 
then it is weakly P-contractive. 
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Proof. Assume that p : X k — > X is a coordinatewise p-contractive /c-mean. We equip 
X k+1 with the sup metric 

d((x 1: . . . , x fc+ i), (yi, . . . , j/fc+i)) := max{rf(x i , y 5 ) : 1 < j < k + 1}. 

We show by induction on n that for any x £ X fc+1 and any two adjacent coordinates 
(/3"(x)) i and (/? n (x)) i+1 , 

d((/3 n (x)),, (/T(x)) i+1 ) < ^(xi.Xi+i). 

For n = 1, we note by the coordinatewise p-contractive property that 

<W( X ))»» (/3(x)) i+i ) = d(//(7r^(x)),//(7r^+i(x)) < 

since vr^j(x) and 7iy i+1 (x) differ only in the ^-coordinate, where they have the en- 
tries Xi + \ and Xi resp. Assume the validity of the inductive hypothesis for n. Since 
/3 n+1 (x) = /3(/3 n (x)), we have from the case n — 1 that 

d((r +1 (x)),, (/T +1 (x)) i+1 ) < p ■ d((/r(x))i, 03 n (x)) m ). 

By the inductive hypothesis, the latter is less than or equal p • p n d(xi, Xj+i) = 
p n+1 <i(a;j, This completes the induction. We then conclude from the tri- 

angle inequality that since any two entries of /3 n (x) are at most k steps apart, 
we have d(p n (x) h p n (x) j ) < kp n A(x), and thus A(/? n (x)) < fcp n A(x). Therefore 
lim n A(/3 n (x)) = 0. □ 

Note that if is power convergent, then p must be weakly /3-contractive. The 
next proposition provides a converse. 

Proposition 3.9. Let X be a complete metric space endowed with a weakly j3- 
contractive k-mean p. If X is uniformly locally convex, then (3 is power convergent, 
so that a f3 '-extension exists. 

Proof. For x £ X, set C n (x) equal to the closed convex hull of |/9 n (x)|. By hypothesis 
A(/3 n (x)) = diam|/? n (x)| — > and then by uniform local convexity diam C n (x) — > 0. 
Note that since C n (x) is convex, it contains |/? m (x)| for all m > n, and hence contains 
C m (x). Thus the collection {C n (x)} is a decreasing sequence of closed convex sets 
whose diameters converge to 0. Since X is a complete metric space the intersection 
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consists of a single point {x*}, and it is now easy to show that /3 n (x) converges to 
the (k + l)-tuple with all entries x*. □ 

We next single out another property that will be useful for inductively /^-extending 
a mean and give some useful equivalences. 

Lemma 3.10. Let X be a metric space endowed with a k-mean ji. We endow X k 
and X k+1 with the sup metric that takes the supremum of the distances between each 
of the corresponding coordinates. Then the following three conditions are equivalent: 

(1) for all x = (x 1: . . . , x k ), y = (yi, . . . , y k ) G X k , 

d(fi(x),(j,(y)) < m.a^{d{x j ,y j ) : 1 < j < k}; 

(2) The mean /i : X k — > X is Lipschitz with Lipschitz constant 1 {hence, in 
particular, is continuous) . 

(3) The map j3 : X k+1 — > X k+1 is Lipschitz with Lipschitz constant 1. 
These conditions imply 

(4) X is closed ball convex. 

Proof. (1) (2): the right-hand side of (1) is the definition of the sup metric, so the 
two statements are equivalent. 

(2) =>- (3): In each coordinate the map (3 is a projection followed by fi, a composition 
of maps with Lipschitz constant 1, and thus has Lipschitz constant 1. Since this holds 
in each coordinate, it holds in the sup metric. 

(3) (2): (3) (2): Fixing some z G X, we have for x G X k that /i(x) = 
■Ki(/3(z, x)), and the right-hand side is a composition of maps of Lipschitz constant 1. 

(2) =>• (4): For e > and x G X, y u . . . , y k G X k , we have for y := (y u ...,y k ) 

d(x,fx(y)) < d(n(x, . . . ,x),n(y)) < d((x, ...,x),y) = maxd(x,yi) < e 
provided d(x, yi) < e for all i. Thus X is closed ball convex. □ 

Definition 3.11. A /c-mean ji on a metric space X is called nonexpansive if it satisfies 
for all x = (xi, . . . , x k ), y = (y u . . . , y k ) G X k , 

(3.1) d(n(x),(i(y)) < max{rf(x i ,i/ i ) : 1 < j < k}, 

or equivalently condition (2) or (3) of the preceding lemma. 
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Lemma 3.12. If p is a nonexpansive k-mean on a metric space X and if p has a 

(3 '-extension p, then fi is nonexpansive. 

Proof. Let tci : X k+1 — > X denote projection into the first coordinate. For x e X k+1 , 

/2(x) = 7Ti(Hm/? n (x)) = lim(7ri o/T)(x). 

n n 

Since /i is the pointwise limit of Lipschitz maps 7Ti o /5 n of Lipschitz constant one (by 
Lemma [3. 101) . it is also. □ 

The next proposition is the principal tool that allows us to extend means inductively 
to higher order. 

Proposition 3.13. Let X be a complete metric space equipped with a nonexpansive, 
coordinatewise p-contractive (0 < p < 1) k-mean p : X k — »■ X, k > 2. Then 
the bary centric operator (3 is power convergent, and hence there exists a (unique) 
continuous (k+l)-mean p : X k+1 — > X that (3-extends p. Furthermore, p : X k+1 — > X 
nonexpansive and coordinatewise p-contractive. 

Proof. By Lemma 13.81 /i is weakly contractive. Since p is nonexpansive, by Lemma 
13.101 X is closed ball convex, hence uniformly locally convex (Lemma 13.41) . and thus 
(3 is power convergent and has a /5-extension to a (k + l)-mean /i by Proposition 13.91 
By Lemma [3. 121 the mean p is nonexpansive and hence continuous (Lemma 13.10( 2)). 

To finish we show that p is coordinatewise p-contractive. Let x, y G X k+1 differ 
only in the j^-coordinate, Xj ^ jjj. Then by coordinatewise p-contractivity 

d{(P(*))i> (P(y))i) = d W^( x )).^(^(y))) <pd(xj,yj), 

since vr^j(x) and vr^j(y) differ in at most one coordinate, and are then Xj and yj in that 
coordinate. Since the inequality holds for each i, we have d(/3(x), /3(y)) < pd(xj,yf). 
Since /3 is nonexpansive by Lemma [3. 101 we conclude that <i(/3 n (x), /3 n (y)) < pd{xj,yf) 
for all n. Taking the limit as n —>■ oo, we obtain d(p(x) , p(y)) < pd(xj,yj). □ 

The next theorem is the culmination of this section. It follows from a straightfor- 
ward induction using the preceding proposition. 

Theorem 3.14. Let X be a complete metric space equipped with a nonexpansive, 
coordinatewise p-contractive (0 < p < 1) k-mean p : X k — >• X , k > 2. Then there 
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exists uniquely a family of continuous means p n : X n — > X, one for every n > k, such 
that each is a ^-extension of the previous one. Furthermore, each p n is nonexpansive 
and coordinatewise p- contractive. 

Example 3.15. Consider on R the mean p(x,y) = sx + (1 — s)y, where < s < 1. 
Set p = max{s, 1 — s}. Then it is an elementary calculation to verify that p is 
coordinatewise p-contractive and nonexpansive. Hence p inductively /3-extends to an 
n-mean for all n > 2. For example, if m(x, y) = (2/3)x + (l/3)y, then one varifies that 
m(x,y,z) = (2/5)a; + (7/20)y + (1/4)2 is a /5-invariant extension of m, and hence must 
be its three- variable /^-extension. If one uses the alternative barycentric operator (3*, 
then one obtains the extension m 3 (x,y,z) = (4/7)x + (2/7)y+ (l/7)z. 

4. Convex means 

In general a metric space may have none, one, or many midpoints between two given 
points in the space. (Recall that m is a midpoint of a and b if d(m, a) = d(m, b) = 
(l/2)d(a, b).) We wish to consider the setting where possibly many midpoints may 
exist, but there is a distinguished midpoint, and these distinguished midpoints appear 
in a "convex" manner. 

Definition 4.1. A symmetric mean p : X x X — > X, written p(x,y) = x#y, on 
a complete metric space X is called a convex mean if it satisfies the basic convexity 
condition 

(4.1) d(x#z, y#z) < ]-d(x, y) for all x,y,z e X. 

Lemma 4.2. For a convex mean, x#y is a metric midpoint for allx,y. 

Proof. By the basic convexity condition d(x#y, y = yj^y) < (l/2)d(x, y) and similarly 
d(x#y,x) < (l/2)d(x,y). Thus 

d(x, y) < d(x, x#y) + d(x#y, y) < ^d(x, y) + ^d(x, y) = d(x, y). 

It follows that each inequality is an equality, so d(x, x#y)+d(x#y, y) = d(x, y). Hence 
adding together the two inequalities in the first line of the proof gives an equality, so 
each inequality is an equality. □ 
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Note that in the case there is only one metric midpoint between two points x, y it 
must be that case that x#y is that midpoint. 

The next proposition gives a useful equivalence for convexity. 

Proposition 4.3. Let (X,d) be a complete metric space equipped with a mean \i. 
Then \x is a convex mean if and only if 

d(x#y, u#v) < -d(x, u) + -d(y, v) for all x, y,u,v G X. 

Proof. For a convex mean d(x#y, u#y) < (l/2)d(x, u) and d(x#y, u#v) < (l/2)d(y,v). 
The condition of the theorem now follows by adding the equations and an application 
of the triangle inequality. 

The reverse implication follows by choosing, x = x, z = y = v, and y = u. □ 

Proposition 4.4. A convex mean inductively j3-extends to a symmetric, nonexpan- 
sive, coordinatewise (1/2) -contractive n-mean for every n > 2. 

Proof. Note that the definition of a convex mean is that of a symmetric coordinatewise 
(l/2)-contractive 2-mean. Proposition 14.31 further yields that it is nonexpansive, since 

d(u#y,u#v) < -d(x,u) + -d(y,v) < m&x{d(x,u),d(y,v)}. 

Thus by Theorem 13. 141 we obtain inductively a /^-extension for every n that is nonex- 
pansive and coordinatewise (l/2)-contractive. By Proposition 12.41 each extension is 
symmetric. □ 

Example 4.5. Let X be a Banach space (or a closed convex subset thereof) and define 
the symmetric 2-mean fi(x,y) = (l/2)(x + y). This is the midpoint with respect to 
the norm metric, and is easily seen to be a convex mean. Setting /ifc(xi, . . . , x^) = 
=1 Xi, one verifies directly that fik+i is the /3-extension of fik, so /i inductively 
/3-extends to the standard arithmetic mean fi n for all n. 

5. Hadamard spaces 

A metric space X is said to satisfy the semiparallelogram law if for any two points 
Xi,X2, G X, there exists z6l that satisfies for all 

d(x\, x 2 ) 2 + 4:d(x, z) 2 < 2d(x, Xi) 2 + 2d(x, a; 2 ) 2 . 
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It follows readily that z is the unique midpoint between x\ and x^. A Hadamard space 
(occasionally called a Bruhat-Tits space) is a complete metric space that satisfies the 
semiparallelogram law. 

Using a metric notion for an upper bound of curvature (geodesic triangles in the 
metric space satisfy certain inequalities when compared with test triangles), one calls 
a metric space a CAT(K,)-space if it is a geodesic space (each pair of points can be 
connected by a metric geodesic) satisfying the curvature bound condition for the 
real number k (see (2J Chapter i] or [51 Chapter II. 1]). The CAT(0)-spaces are the 
non-positively curved spaces. A metric space has an alternative characterization as a 
Hadamard space: it is a simply connected, complete, geodesic CAT(0)-space (see [21 
Proposition 5.1, Chapter 1] or [5, Exercise 1.9, Chapter II. 1]). 

What is important for our current purposes is the following remark: 

Remark 5.1. Let A be a Hadamard metric space and define a 2-mean by n(x,y) is 
the unique midpoint between x and y. This defines a convex mean in the sense of the 
preceding section, see [21 Proposition 5.4, Chapter I] or [51 Proposition 2.2, Chapter 
II. 2]. Hence by the preceding section this mean may be /3-extended to an n-mean for 
every n > 2. 

A wide variety of Hadamard spaces and constructions for new Hadamard spaces 
from old appear in [2] and [5]. Some examples include Hadamard manifolds (simply 
connected complete Riemannian manifolds with nonpositive sectional curvature), par- 
ticularly simply connected symmetric spaces of noncompact type, finite-dimensional 
hyperbolic geometries over the reals, complexes, and quaternions, symmetric cones, 
Tits buildings, and various examples obtained by coning and gluing. 

Of particular interest to us is the example of the manifold of positive definite 
matrices endowed with the usual Riemannian metric called the trace metric. This 
metric yields a Hadamard manifold and the midpoint mean operation in this case is 
precisely the geometric mean of the two positive definite matrices; see [12] and the 
references there. Using the fact that the length metric satisfies the semiparallelogram 
law, hence is a convex metric with the midpoint operation being a convex mean, we 
obtain the following alternative derivation of the principal result of [I]: 
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Corollary 5.2. Let X denote the set of positive definite real or complex matrices 
equipped with the Riemannian trace metric. Then the midpoint operation for the 
corresponding length metric, which is precisely the geometric mean, defines a convex 
2-mean, which (by Proposition 4-4) (3-extends to an n-mean for each n > 2. 



6. Iterated means 

A standard construction technique for means is iteration, the arithmetic-geometric 
mean being the best known example. In this section we develop machinery for showing 
that certain iterated means are coordinatewise p-contractive and nonexpansive, hence 
admit /3-extensions of all orders. We apply this machinery to operator means in a 
later section. 

Definition 6.1. Let A, v be 2-means on a complete metric space X. Starting with 
Ai = A and v\ = u, we define inductively the sequences of means {A n } and {z/ n } by 

K+i(x,y) = X(X n (x,y),u n (x,y)), u n+1 (x,y) = v(\ n (x, y),u n (x,y) ) 
(resp. u n+1 (x,y) = u(X n (x, y), u n (x, y)), X n+1 (x,y) = X(X n (x,y), u n+1 (x,y)). 

If there exists a 2-mean \x such that for all x, y G X, lim n X n (x, y) = /j,(x,y) = 
]im n i> n (x, y), then \i is called the iterated composition (resp. skewed iterated compo- 
sition) of A and v and denoted fj, — A * v (resp. /x = A * s v). 

We begin with a useful lemma that ensures convergence. 

Lemma 6.2. Let {x n }, {y n } be sequences in a complete metric space X satisfying 
one of the two following conditions: 

(i) for each k > 1, xt+i is a midpoint of Xk andyu and d(xk+i,yk+i) < d(xk+x,yk); 
or 

(ii) for each k > 1, Xk+i is a midpoint of Xk and yu+i and d(xk,yk+i) < d(xk,yk)- 
Then both sequences are Cauchy and converge to the same point. 

Proof. Assume (i). For any n > 1, we have by hypothesis 

d(x n+ x,y n+1 ) < d(x n+ x,y n ) = ^d(x n ,y n ), 
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where the last equality follows from that the fact that x n+ i is a midpoint of x n ,y n . 
Similarly 

It follows by induction resp. induction and the triangle inequality that 

1 1 \ 1 

d(x n ,y n ) < -^— [ d{x u y l ) resp. d(x n+k , x n ) < I — ) d ( Xu y ^ < ^i^ 1 ' 

Thus the sequence {x n } is Cauchy, and hence converges, and the sequence {y n } must 
also approach the same limit. 

Part (ii) follows by applying part (i) to the sequences {x n } and {z n }, where z n = 

Vn+l- 

□ 

Proposition 6.3. Let X be a convex mean and v be a nonexpansive mean on a 
complete metric space X . Then the iterated composition fi = A * v resp. the skewed 
iterated composition fi = A * s v exists and is nonexpansive. 



Proof. For x, y e X, we set x n = X n (x, y) and y n = u n (x, y) (see Definition l6.ll) . Then 
x n +\ = K x n,y n ) and y n+ i = u(x n ,y n ). We observe that 

d(x n+ i,y n+ i) = d(i/(x n+ x,x n+ i),v(x n ,y n )) < max ),d(x n+1 ,y n )}, 

where the last inequality follows from the fact that v is nonexpansive. Since A is 
a convex mean x n+ \ = \(x n ,y n ) is a midpoint for x n and y n , hence d(x n+ \, x n ) = 
d(x n+ i,y n ), and thus d(x n+ i, y n +i) < d(x n+ i,y n ), i.e., condition (i) of Lemma [6T21 is 
satisfied. It thus follows that lim„ x n = lim n y n exists, and we define this limit to be 
fi(x,y). If x = y, then it is immediate that x = x n = y n for all n, so fi is a mean. 
Thus the iterated composition = A * v exists. 

For the case of the skewed iterated mean, we set X\ = A (2, y), yi— v{x, y) and 

Z/n+i = v(x n ,y n ), x n +i = \(x n ,y n +i). 

Then 



d(x k ,y k+1 ) = d(v(x k ,x k ),v(x k ,y k )) < m&x{d(x k ,x k ),d(x k ,y k )} = d(x k ,y k ), 
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where the inequality follows from the nonexpansive property. Thus Lemma 16.2( h) is 
satisfied. That the skewed iterated composition p = X * s v exists now follows as in 
the preceding paragraph. 

It follows from Proposition 14. 31 that A is Lipschitz with Lipschitz constant 1 and the 
same holds for v since it is non-expansive. In both cases the higher numbered means 
X n and v n are built up from these by products and compositions, so are also 1-Lipschitz 
(recall that product metrics are always the sup metric). Since p is the pointwise limit 
of the sequence {A n } (and {v n }), it is also 1-Lipschitz, i.e., nonexpansive. □ 

Proposition 6.4. Suppose in a complete metric space X that X is a convex mean 
and v is coordinatewise p' -contractive, < p' < 1, and nonexpansive. Then the 
iterated composition X * v resp. the skewed iterated composition X * s v exists and is 
coordinatewise p- contractive, p = max{l/2, p'}, and nonexpansive, hence (3-extends 
to all orders greater than two. 

Proof. By Proposition 16.31 the iterated composition and skewed iterated composition 
both exist and are nonexpansive. 

We establish that p = X * v is coordinatewise p-contractive. Let a, b, c G X. To 
calculate p(a, b) and p(a, c) we define inductively 



b x = X{a, b), bf = u(a, b), b k+1 = X(b k , &+), = u{b k , &+) 
Ci = A(a,c), cf = is(a,c), c k+1 = A(c^,c+), c+ +1 = u(c k ,c+). 



Note that b k = Xk(a,b), b k = z/fc(a,6), c k = Afc(a, c), c k = z/fc(a, c). We have 
d(bi , Ci) = d(X(a, b), A(a, c)) < (l/2)<i(6, c) by convexity of A and similarly d(bf, c^) < 
pd(b, c) by coordinatewise p-contractivity of v. 
We claim that by induction 



d(b n ,c n ) <pd(b,c) and d(&+,c+) <pd(b,c). 
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By the preceding paragraph it holds for n — 1. Assume that it is true for n = k. 
Then 

d ( b k+n c k+i) = rf (A(^,^),A(Cfc,c+)) 

< d{X{b k , bt), X(c k , bt)) + d(X(c- k , bt), X(c- k , c + k )) 

< \d{b- k ,c- k ) + l -d{bt,ct) 

< -(pd(b,c) + pd(b,c)) = pd(b,c). 
Using the nonexpansivity of u k+ i, we obtain 

d ( b t+v4 + i) = d n b k^ b t)^( c k^ c t)) 

< max{d(b k ,c k ),d(bt,ct)} 

< max{prf(6, c), pd(b, c)} = pd{b, c). 

This completes the induction. Note that in the alternative notation we have shown 
that d(X n (a, b), X n (a, c)) < pd(b, c) and d{y n {a, b), v n {a, c)) < pd{b, c) for all n e N. 

Since p — X * v, lim n b~ = lim n A n (a, b) = p(a, b), lim n c~ = lim n A„(a, c) = p(a, c). 
By continuity of d(-, ■) and the preceding paragraph, it follows that d(p(a, b), p(a, c))) < 
pd(b, c). 

The proof that the skewed iterated composition p = X * s v is coordinatewise p- 
contractive is similar, but contains a twist or two. To calculate p{a, b) and p{a, c) for 
a,b,c G X, we define inductively 

b+ = u(a,b), b± = X(a,b), bt +1 = v(b k ,bt), b k+1 = X(b k ,bt +1 ) 
c + = u(a,c), c{ = X(a,c), ct +1 = u{c k ,4), c k+1 = X(c k ,ct +1 ). 

We have d(b^,Ci) = d(X(a,b), X(a, c)) < (l/2)d(b,c) by convexity of A and similarly 
d{pt, cf) < pd{b,c) by coordinatewise p-contr activity of v. 
We claim that by induction 

d(b~, c~) < pd(b, c) and c+) < pd(b, c). 
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By the preceding paragraph it holds for n = 1. Assume that it is true for n = k. 
Using the nonexpansivity of u, we obtain 

< max.{d{bj~,cj;) } d(b£,c£)} 

< max{pd(b, c),pd(b, c)} = pd(b, c). 

It then follows that 

d ( h k+n c k+i) = d(\(b k ,b+ +1 ),\(c k ,4 +1 )) 

< d(\(b k , b+ +1 ), \(c k , b+ +1 )) + d(\(c-, 6+ +1 ), A(cjr, c+ +1 )) 

< ^( & fc> c fe) + ^( & fc+i> c fc+i) 

< ^(pd(b,c) + pd(b,c)) = pd(b,c). 

This completes the induction. 

By hypothesis lim n 6~ = lim n A Tl (a, b) = p(a,b), lim n c~ = lim n A n (a, c) = p(a,c). 
By continuity of d(-, •) and the preceding paragraph, it follows that d(p(a, b), p(a, c))) < 
pd(b, c). 

The last assertion of the the proposition now follows from Theorem 13.141 

□ 

7. Categorical constructions 

In this section we consider the behavior of mean extensions with respect to standard 
constructions such as continuous images, products, and subspaces. 

Definition 7.1. A function g : (X, p) — » (Y,v), where p,u are fc-means on X and 
Y respectively, is called a k-mean homomorphism or homomorphism for short, if 
g o p = v o g k) that is, the following diagram commutes 

X -U Y 



X k 9U > y k 

where g k : X k -> F fc , ...,x k ) := (#(xi), . . . , #(x fc )). 
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Proposition 7.2. Let {X,jx) and (Y,u) be topological k-means, and let g : X — > Y 

be a continuous k-mean homomorphism. 

(i) If each of \i and v (3-extend to (k + l)-means fx and v resp., then g : X — ► Y 
is a (k + l)-mean homomorphism. 

(ii) If g is surjective and fi /3-extends to a (k + l)-mean fx, then v /3-extends to a 
(k + l)-mean v, and g is then a (k + \)-mean homomorphism. 

Proof, (i) It follows directly from the fact that g is £;-mean homomorphism that 
PyQu+i = 9k+iPx '■ X k+1 — > Y k+1 (indeed commutation of with f3 is an equiva- 
lence). By induction 

ft 9k+1 = g k+ i(3 n x ; x k+l - Y k+1 

for all n > 0. Taking the limit of both sides as n — ► oo and projecting into the first 
coordinate yields (i). 

(ii) For y £ Y" fc+1 , there exists x £ X k+1 such that ^ + i(x) = y. Again we have 

0?(y) = ^f/ fc +i(x) = g k+ iP n x (x) 

for n > 0. By hypothesis the right-hand side converges to a diagonal element with 
entries g(Jl{x)) as n — > oo, so that the left-hand side also converges to a diagonal 
element. Thus i/ /3-extends to P. The last assertion follows from (i). □ 

Example 7.3. Let 

be a quasi- arithmetic mean defined on IR + by a continuous strictly monotonic function 
/ [IT] - By the preceding proposition applied to g — f^ 1 and Example 14.51 



1 k 

l!l,{Xl,...,X k ) = f 1 (T^2f{Xi))' 



i=l 



Note that the arithmetic, geometric, and harmonic means belong to this class by 
taking the identity map, the logarithmic map, and the inversion map (on the positive 
reals) respectively. More generally one can take on the positive reals the generalized 



or power mean m(x,y) = ((x a + y a )2) 1 ^ a for with f{x 



x 
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An analogous construction and characterization of the higher order means remains 
valid for the power means on the space of positive definite matrices. Note that the 
case a = 1 gives the arithmetic mean and the case a — — 1 gives the harmonic mean. 

Definition 7.4. Let p, v be /c-means on X and Y, respectively. Define p x v : 
(X x Y) k -> X x y by 

(// x u)(x 1 ,y 1 ,x 2 ,y 2 ,...,x k ,y k ) = (p(x 1: x 2 , . . . ,x k ),v(y 1: y 2 , . . . ,y k )). 

Then /x x z/ is a /c-mean on X x Y, called the product mean of p and za Indeed, if 
z = (x,y,x,y, . . . ,x,y) G (X x y) fe for fixed G X x y then (// x i/)(z) = 

(fi(x, x,...,x), u(y, y,...,y)) = (x, y). 

Theorem 7.5. Let (X, d\) and (Y, d 2 ) be complete metric spaces equipped with non- 
expansive, coordinatewise p-contractive resp. p' -contractive (0 < p, p' < 1) k-means 
p : X k — > X resp. z/ : y fe — > y. T/ien i/ie mean p x u is a nonexpansive, coordinate- 
wise max{p, p'} -contractive k-mean on X x Y. Furthermore, for n > k, its n-mean 
[3 -extension (p x u) n coincides with the product mean p n x u n of the individual f3- 
extension n-means: 

(p x v) n = p n x v n . 

Proof. For z = (x 1 ,y 1 ,x 2 ,y 2 , . . .,x k ,y k ) G (XxY) k , we denote z x = (x 1 ,x 2 , ...,x k ) G 
X k , z y = (y 1 ,y 2 ,...,y k ) G Y k . Then (p x i/)(z) = (p(z x ),u(z y )). 

Let z = (x n ,y n ,x 12 ,y 12 , . . . ,x lk ,y lk ), w = (x 21 , y 21 , x 22 , y 22 , . . . , x 2k , y 2k ) G (X x 
Y) k . Then by nonexpansive property of p and v, 

d((px v)(z),(px = d((p(z x ),is(z y )),(p(w x ),v(Wy)j^ 

= max{d 1 (p(z x ),p(w x )),d 2 (iy(z y ),iy(w y ))} 

< ma^{maix{d 1 (x lj ,x 2j )},maix{d 2 (y 1 j,y 2j )} : 1 < j < k} 

= max{d 1 (x lj ,X2j),d2(yij,y2j) ■ 1 < j < k} 
= max{d((a:i j , y Xj ), (x 2j ,y 2j )) : 1 < j < k}, 

which implies that p x v is a non-expansive /c-mean on X x Y equipped with the 
sup-metric. 
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If z and w differ only in one coordinate of (X x Y) k , say (x±j, y±j) ^ (x 2 j, y 2 j), but 
(xii,yu) = (x2i,V2i), 1 < i 7^ j < k, then the inequality in preceding argument turns 
into 



d^fj, x v)(z), (p x p)(w)J = d[[p J {z x ),v{z y )), (fJ,(w x ), u(w y )) 

= m&x{d 1 (ii(z x ),fi(w x )),d 2 (v(z y ),v( y w y ))} 

< max{pdi(xi i; x 2j ), p'd 2 (yij, y 2 j)} 

< max{p, p'} max{d 1 (x lj ,x 2j ), d 2 (y lj7 y 2j )} 

< max{p,p }d{{x lj ,y lj ), {x 2j ,y 2j )j. 

Therefore p x v is a coordinatewise max{p, p'}- contractive fc-mean on X x K 

Next, we will prove (p x z/)„(z) = (p n (z x ), v(zy)), z G (X x F) n , n > A; by induction. 
The case n = k follows by the definition of p X v. Suppose that the assertion holds 
true for n — 1. Let 7 : (X x Y) n — > X x K be defined by 7(z) = (p n (z x ), v n (z y )). 
Then 7 is continuous and hence it suffices to show 7 o /3 n = 7 by the uniqueness of 
mean extension (Proposition 12. 4p where /3„ is the barycentric operator on (X x Y) n 
obtained from the (n — l)-mean (p x v) n -\. 

For z = (zi, z 2 , . . . , z n ) G (X x Y) n , Zi = (xj, yi) G X x Y, we have from induction 

and then from /3 n (z) = ((p x z/^^iz), (p x z/) n _ 1 (7T^ 2 z), . . . , (p x z/) n _i(7iy n z)Y 

(3 n (z) x = ^p n -l(7r^lZ x ),p n _i(7r^2Za : ), • • • ,Pn-l(71>nZx.)) = Ax( z z) 
AjX 2 )^ = ^n-l(^lZ 2/ ),^n-l(7r^2Z 2/ ),...,^ rl -l(7r^„Z 2/ )J = & 



'vV'y) 



where P^,P U denote the barycentric operators of the n- means p n , v n \ p n o (3^ = p n 
and v n o (3 U = v n . Now, 

(7o/3 n )(z) = 7(/3 n (z)) = (pn(/5r i .(z) x ),z/ n (/5 n (z) ?/ )) 

= (A i n(/5 M (z x )),z/ n (/3 1 ,(zj / ))j 
= (/x n (z iE ),i/ n (Zj / )) = 7(z) 
which completes the claim. □ 
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Remark 7.6. The product mean satisfies the associative law 

H x (v X u) = (/! x v) X UJ 

for any fc-means /i, u, and uj on X, Y, Z. If these are nonexpansive and coordinatewise 
contractive, then its n-mean extension satisfies 

(/i X V X kj) n = fi n X Z/„ X 

The next result of the section is quite straightforward and hence the proof is omit- 
ted. 

Proposition 7.7. If Z is a nonempty closed k-submean of a topological k -mean (X,n) 
that ^-extends, then Z also (3-extends and fi\Z k+1 = 

Corollary 7.8. Suppose that a k-mean (X, /i) restricts to a nonexpansive, coordi- 
natewise p n - contractive mean on A n , where < p n < 1 for each n and A n is an 
increasing sequence of closed convex sets with X = [j n A n . Then fi inductively f3- 
extends to j-mean for each j > k in such a way that the restriction to A n is the 
appropriate f3-extension of the restriction of fi to A n . 

Proof. Given any j-tuple in X 7 for j > k, there exists some A m that contains all entries 
of the tuple. Applying Theorem 13.141 to the restriction of fi to A m , we conclude that 
A*U m /3-extends for every index greater than k. Thus in particular the appropriate 
extension exists to evaluate the given j-tuple. It is clear that if a larger A m+n is 
chosen, one obtains the same calculaton. Thus the /3-extension is independent of the 
containing A m , and hence we obtain a /3-extension of /i on all of X. □ 

8. Stability and reductions 

Definition 8.1. A k-mean pna set X is called ^-stable if the graph of \i is invariant 
under j3, that is, 

VX G X fc+1 ,7T fe+ i(x) = ^(7T^ fc+ i(x)) 7T fc+ iG0(x)) = /I (/? (x) ) ) . 

A (k + l)-mean v on X is called a stable extension of fi if 

u(x!, ...,x k , n(xi, x k )) = . . . ,x k ) for all x u . . . , x k e X. 
Conversely \i is called a stable reduction of v. 
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In [9] Horwitz says that v is type 2 invariant with respect to p if v is a stable 
extension of p. He only considers the case of a 2-mean p and a 3-mean v. 

Proposition 8.2. // a topological k-mean p is (3-stable and admits a (3-extension p, 
then p is a stable extension of p. 

Proof. Note that (3 does not change the last coordinate of any k + 1-tuple x = 
(xx, . . . . . . ,Xk)) in the graph of p. Hence if p is /3-stable, it follows that 

/3™(x) has the same last coordinate for all n. Since (3 power converges to p, it fol- 
lows that the diagonal limiting value has entries the last coordinate of x, namely 
. . . , Xk). Thus /i(x) = fi(xi, . . . , Xk), so p is a stable extension of p. □ 

Proposition 8.3. If X is a complete metric space equipped with a coordinatewise 
p-contractive (symmetric) {k + l)-mean v , k > 2 then v admits exactly one stable 
reduction (which is also symmetric). 

Proof. Suppose that X is a complete metric space equipped with a coordinatewise p- 
contractive (k + l)-mean v. The map g : X — > X defined by g(x) = u(xi, . . . , x^, x) is 
by hypothesis p-contractive, and hence has a unique fixed point. Define p(xi, . . . , Xk) 
to be this fixed point. It follows immediately that p is a stable reduction of u, and 
uniqueness of the fixed point guarantees that it is the unique reduction. Note that 
the fact that v is a mean and the definition of p imply that p is also a mean. The 
symmetry of p follows directly from that of v. □ 

The following corollary is an immediate consequence of Propositions I3.13[ I8.2[ and 

Corollary 8.4. Let X be a complete metric space equipped with a nonexpansive, 
coordinatewise p-contractive k-mean p. If p is [3-stable, then the unique f3-extension 
p is a stable extension of p, and p is the unique stable reduction of p. 

Definition 8.5. A 2-mean p on a set X satisfies the limited medial property if 
p(a,b) = p(x,y) =: m implies that p(p(a,x), p(b,y)) = m. For m E X, 

X m := {(x,y) E X 2 : p(x,y) = m}. 
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Remark 8.6. (1) If a 2-mean fi(x,y) = xj^y satisfies the limited medial property, 
then it is /3-stable since 

p(x,y,x#y) = (y#(x#y),x#(x#?/),x#y). 

The latter is in the graph of /x, since by limited mediality x#y = (x#y)#(x#y) 
implies (x#(x#y))#(y#(x#y)) = x#y. 

(2) It was shown in [lj that the matrix geometric mean for positive definite matrices 
satisfies the limited medial property. This was extended to very general notions of 
geometric mean in [15], in particular for the geometric mean of positive operators on 
a Hilbert space or, more generally, for the positive elements of a C*-algebra. Hence by 
part (1) and the earlier results each fc-extension of the geometric mean yields both the 
higher order ones by stable extension and the lower order ones by stable reduction. 

Lemma 8.7. The mean \x satisfies the limited medial property if and only if (/x x 
/x)(X m x X m ) C X m for each m G X. 

Proof. Let (a,b), (x,y) G X m , that is //(a, b) = fi(x,y) = m. Then (/x(a, x), /x(o, y)) G 
X m (limited medial property) if and only if/x(/x(a, x),/i(b, y)) = (/xx/x)((a, b), (x, y)) = 
m. □ 

Proposition 8.8. Let X be complete metric space equipped with nonexpansive, coor- 
dinatewise p-contractive 2-mean fi(x,y) = xj^y satisfying the limited medial property. 
Ifxi#Vi = m for all I <i < n, then /i n (x)#/i n (y) = m, where x = (x t , . . . ,x n ),y = 
(yi,...,y n ) eX n . 



Proof. By previous lemma (/x x /i)(X m x X m ) C X m , and X m is a non-empty closed 
subset of X 2 . Let uj := (/i : 
X^ -> X m of lj is given by 



subset of X 2 . Let lj := (/i x //) . By Theorem 17.5( 3). the n-mean extension lj h : 



Suppose that x^yi = m, i = 1, 2, . . . , n. Then z := (x u y u x 2 , y 2 , • • • , x n , y n ) G X™ 
and so u n (z) = (/i n (x), /x n (y)) G X m which implies that /i n (x)#/i n (y) = m. □ 



Remark 8.9. Under the assumption of the preceding proposition, if aj^-b = xj^y = m 
then X m C {(2, w) G X x X : 7x3(0, 6, z)#fi 3 (x, y, w) = m}. 
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9. Ordered convex metric spaces 

Throughout this section we assume that X is a complete metric space equipped 
with a non-expansive, coordinatewise p-contractive /c-mean p : X k X,k > 2, and 
p n (n > k) denotes the non-expansive, coordinatewise p-contractive n-mean obtained 
inductively. We further assume that X is equipped a closed partial order < . We let 
<k be the product order on X k defined by 

(xi, ...,x k ) <k (yi, ■■■,yk) if and only if Xi < y h 1 < i < k. 
Recall for a map g : X — > Y, we let g k : X k — > Y k , (x 1: . . . , x k ) i— > (g(xi), . . . , g(x k )). 

Definition 9.1. A /c-mean z/ on X is said to be monotone for the partial order order 
< if z/(x) < z/(y) for any x, y G X k with x < y. 

Theorem 9.2. (1) If a nonexpansive, coordinatewise p-contractive k-mean p is mono- 
tone for the closed partial order <, i/ien p n is also monotone for any n > k: 

X < n y => Pn(x) < Pn(y)- 

(%) Lei (X, <5) 6e another complete metric space equipped with a nonexpansive, coor- 
dinatewise p-contractive k-mean v : X fe — > X. Suppose that < is closed in the product 
topology (X, d) x (X, 5) and either p or v is monotone with respect to < . Then p < v 
implies p n < v n for any n > k. 

Proof. (1) We induct on m beginning at k. Suppose that p m is monotone. Let x, y e 
X m+1 with x < m+ i y. Then 7r^x < m 7r^y for all 1 < i < m + 1. By the induction 
hypothesis, p m (vr^x) < p m (ir^y) for each i, and hence /3 m+ i(x) < m+ i /3 m+ i(y). By 
repeated application, we have /5™ +1 (x) < m+ i /3™ +1 (y), n = 1, 2, .... By the closeness 
of the order, lim^oo Pn +1 (x) < n+ i lim^oo /3£ +1 (y) . In particular, p n+ i(x) < p„+i(y). 

(2) Assume that p m (x) < z/ m (x) for all x e X m for some m > k. Let x e X m+1 . 
Then p m (7r^x) < z/ m (7r^jx) for 1 < j < m + 1 and hence 

(9.1) /3(x) < m+1 a(x) 

where /3 and a are the barycentric operators with respect to the means p and v, 
respectively. We will show by induction that /3 n (x) < m+ i « n (x) for all n. Suppose 
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that it holds true for n. Then 7r^y/3 n (x) < 7r^jO;' 1 (x) for all 1 < j < m + 1. If v is 
monotone, J/ m is monotone by (1) and 

(9.2) z/ m (7r^/3 n (x)) < z/ m (7r^a n (x)), 1 < j < m+ 1. 
Therefore, 

/T +1 (x) = f^+i aC9»(x)) 

(^ m (vr^i/9 n (x)), . . . , z/ m (vr^ m+1 /5 n (x))) 

< m+i (^m(vr^ia"(x)), . . . , z/ m (7r^ fc+ ia; n (x))j 

a(a n (x)) = « n+1 (x). 

In the case that p is monotone, (19.21) changes to 

(9.3) /i m (vr^/? n (x)) < /i m (7r^a n (x)), 1 < j < fc + 1, 
and /3 n+1 (x) < m+i a n+1 (x) from 

/T +1 (x) = /3(/T(x)) 

(/^(tt^/T (x)), . . . , /i m (7r^ m+1 /5 n (x))) 

< m+i (/i m (7r^ia n (x)), . . . , /i m (vr^ m+ ia n (x))) 

induction 

< m+1 (^m(^l« n (x)),...,Z/ m (7T^ m+ ia n (x))) 

a(a n (x)) = a n+1 (x). 
Since the order is closed in the product topology, the inequality holds for limits: 

lim /? n (x) < m+1 lim a n (x). 
In particular, /i m+ i(x) < i/ m+ i(x). □ 

Theorem 9.3. Let X and Y be complete metric spaces equipped with nonexpansive, 
coordinatewise p (resp. p' ')- contractive k-means p and u, respectively. Let < be a 
closed partial order on Y. Suppose that the mean v is monotone for the partial order 
< onY. If g : X — > Y is continuous and satisfies g o p < v o g k [resp. g o p > v o g k ), 
then for any n > k, 

g ° Pn < v n ° 9n (resp. g o p n > v n o g n ) . 
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Proof. Suppose that g o // < v o g k ; the other case is similar. Suppose that g o /i m < 
v m ° gm holds true for some m > k. Let a be the barycentric operator with respect to 
the mean v on Y. Then g{ji m {j[^p^)^ < v m g m (j[^pdj for x G X m+1 and 1 < j < m + 1, 
and thus 

g> m+ i(/3(x)) = (#(// m (7iyix)), . . . , g-(^ m (7T^ m+ iX))) 

(^(^(^IX)), . . . , Z/ m (# m (7iy m+ ix))) 
= (z/ m (7T^i5( m+1 (x)), . . . , ^ m (vT^ m+ i5( m+1 (x))) 

= a(fiWi(x)) 



and therefore g m+ i o (3 < m+ i a o g m+ i. Since z/ is monotone, a is monotone for < m 
(Theorem 19.21) . and therefore 

^ m +i(/3 2 (x)) = (3> m+ i/3)(/3(x)) < m+ i a(5f m+ i(/3(x))) < a 2 (g m+1 (x.)) 

and inductively g m+ i(/?™(x)) < m+ i a n (g m+ i(x))) for all n. Since the order is closed, 

5f m+ i(lim /3 n (x)) < m+ i lim a n (g m+1 (x)). 

n— >oo n^oo 

In particular, g(/i m+1 (x)) < u m+ i{g m+ i (x)). Induction on m yields the theorem. □ 

Corollary 9.4. Lei /x, z/ and uj be nonexpansive, coordinatewise contractive k-means 
on complete metric spaces X, Y, Z respectively. Let < be a closed partial order on 
Z and let g : X x Y — > Z be a continuous function. If g o (// x v) < u o g kl then 
g o (/i x v) n < iu n o g n for any n > k; 

g(li n {xi,x 2 , ■ ■ ■ ,x n ),v n (y 1 ,y 2 , . . . ,j/ n )J < Wnfp(jci,2/i), ^(^2, 2/2), • • .,g(x n ,y n )\. 

Proof. By Theorem 17.51 the product mean fi x i/ is a nonexpansive, coordinatewise 
contractive fc-mean on X x y. By Theorem 19.31 the inequality go [fi x u) < uj o g k can 
be extended to n-means: go(/jx z/)„ < uJ n o g n . From 

(/i X Z/)„(Xi, y U X 2 , 2/2, • • • , ^n, 2/n) = ^n(^l, ^2, • • • , ^n), ^n(j/l, J/2, • • • , 2/n)J 

g n (x 1 ,y 1 ,x 2 ,y 2 ,...,x n ,y n ) = (^(x^yt), g(x 2 ,y 2 ), . . . ,g(x n ,y n )\, 
the proof is completed. □ 
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10. Means on Hilbert space operators 

In this and the next section we apply and extend our preceding results to the special 
setting of positive definite operators on a Hilbert space, in particular to positive 
definite Hermitian matrices (in the case the Hilbert space is finite dimensional). 

For a Hilbert space E, let B(E) denote the set of bounded linear operators, S(E) C 
B(E) the symmetric operators, and Q C S(E) the set of positive definite operators 
on E. We define a closed positive order on S(E) by A < B if B — A is positive 
semidefinite. Note that the identity operator / (and indeed any positive definite 
operator) is an order unit for S{E) (that is, S(E) = {J^ =1 [—nl,nl], where in general 
[A, B] denotes the order interval [A, B] = {X e S{E) : A < X < B}. There is a 
corresponding order unit norm given by 

\\A\\ = mf{t >0: Ae [-tl,tl]}. 

This norm generates the same topology on S(E) as does the operator norm. 
We primarily employ the Thompson (or part) metric on Q given by 

d(A,B) = max{\ogM(A/B),\ogM(B/A)} where M{A/B) := inf{A > : A < \B}. 

A. C. Thompson [18] has shown that Q is a complete metric space with respect to 
this metric and the corresponding metric topology on Q agrees with the relative norm 
topology. We list some additional elementary properties of the Thompson metric. 

Lemma 10.1. The Thompson metric on the set Q of positive definite Hilbert space 
operators satisfies 

(i) d{A + B,A + C) < d{B,C); 

(ii) A 1 < A 2 implies d(A 1 + B,A 1 + C)> d(A 2 + 5, A 2 + C); 

(iii) For r > 0, d(rA,rB) = d(A, B); 

(iv) d(A + B,C + D)< max{d(A, C), d(B, D)}; 

(v) d(A~ l , B~ l ) = d(A, B). 

Proof, (i) There exists r > 1 such that logr = d(B,C). Then B < rC, and thus 
A+B < A+rC < rA+rC = r(A+C), and similarly C <rB implies A+C < r(A+B). 
Hence d(A + B, A + C) < logr = d(B, C). 
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(ii) The proof of (i) remains valid for A > 0, and then (ii) follows by rewriting A2 
as A l + (A 2 - Ax). 

(iii) For r > 0, d(rA,rB) = d(A,B) since scalar multiplication by r is an order- 
isomorphism. 

(iv) Suppose that d(A, C) < d(B, D) = logr. Then B < rD, D < rB, A < rC, C < 
rA, and thus A + B < rC + rD = r(C + D),C + D < rA + rB = r(A + B). Hence 
d(A + B, C + D) < logr = d(B, D). 

(v) This follows from the order reversing property of operator inversion. □ 

The results of Section 6 require forming iterated means from means that are 
on the one hand convex, and on the other coordinatewise p-contractive and non- 
expansive. The convex mean we employ is the geometric mean on Q defined by 
A#B = A 1/2 (A~ 1/2 BA~ 1/2 ) 1/2 A 1/2 (see [12] for a variety of other characterizations). 
It is known that the geometric mean is a convex mean with respect to the Thompson 
metric (see, for example, [7], [2] ). We record this fact. 

Lemma 10.2. The geometric mean Aj^B is a convex mean on Q endowed with the 
respect to the Thompson metric. 

We close this section by establishing that the arithmetic mean and the harmonic 
mean are nonexpansive and coordinatewise p-contractive with respect to the Thomp- 
son metric when we restrict to order intervals. This takes some computation. 

Lemma 10.3. For each n eN, there exists p n , < p n < 1, such that the arithmetic 
mean A(X, Y) = (X + Y)/2 and the harmonic mean H(X, Y) = 2(X~ 1 + Y^ 1 )^ 1 are 
coordinatewise p n - contractive and nonexpansive on the order interval [(l/n)I, nl\ C 
Q, where I is the identity operator. 

Proof. By Lemma Il0.ll (iv), the arithmetic mean is nonexpansive. We seek p = p n 
such that for any A,B,C E [(l/n)J, nl), d(A + B, A + C) < pd(B, C) (note that we 
can drop the factor of 1/2 by Lemma [lO.l( iii)). This is equivalent to 

ma,x{\ogM(A+B/A+C),hgM(A+C/A+B))} < pmax{log M(B/C), log M(C/B))}. 

Note that if B ^ C (the desired inequality is trivially true if B = C), then either 
B ^ C or C ^ B, and then there exists r > 1 such that C < rB, r = M(C/B), and 
logr = d(B, C) (or vice- versa with the roles of B and C interchanged). 
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Suppose now that we can find p, < p < 1, such that for all A,B,C G [(l/n)I, nl], 
we have A + C < r p (A + B), where r = M(C/B) = max{M(B/C), M(C/B)}. Then 

d(A + B,A + C) <\ogr p = p\ogr = pd(B,C), 

which is the coordinatewise p-contractive property. We conclude that to establish 
the coordinatewise p-contractive property, it suffices to show the existence of some p, 
< p < 1, such that for all A, B,C G [{l/n)l,nl], if C < rB for 1 < r < n 2 , then 
A + C < r p (A + B). (Note that we can restrict to r < n 2 since C < nl = n 2 {l/n)I < 
n 2 B.) We establish that this is indeed the case by means of the two following claims. 
Claim 1: For given < p < 1 and r > 1, assume that C + (l/n)I < r p (B + (l/n)I) 
whenever C < rB. Then A + C < r p (A + B) for all A > (l/n)J. Indeed 

A + C = (A--I) + -I + C 

n n 

< (A--I)+r p (-I + B) 

n n 

< r p {A--I)+r p {-I + B) =r p {A + B). 

n n 

Claim 2: There exists p, < p < 1, such that for all 1 < r < n 2 , C + (l/n)I < 
r p (B + (l/n)I) whenever C < rB, B,C G [(l/n)l,nl]. Indeed for any < p < 1, 

C + -I < rB + -I = r p B + (r - r p )B + -I. 
n n n 

For the last two terms we have 

(r - r p )B + -I < (r - r p )(nl) + -I = {(r - r p )n 2 + . 

n n n 

To complete the proof of Claim 2, we need to choose p < 1, but large enough so that 

(r — r p )n 2 + 1 < r p for 1 < r < n 2 . The function f(x) = x p — n 2 (x — x p ) — 1 has 

derivative 

f{x) = px p - 1 - n 2 (l - px"- 1 ) = (1 + n 2 )px p ~ x - n 2 . 

Since the limit of the right-hand expression is 1 as p — > we conclude that the 
derivative is positive for all x G [1, n 2 ] for large enough p below 1. Thus / is increasing 
on [l,n 2 ] for p near, but below, 1, and hence 

(r — r p )n 2 + 1 < r p for any 1 < r < n 2 , \1 — p\ < e 

for some e > 0. 
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The case of the harmonic mean follows from the fact that the inversion is an 
isometry (Lemma 110.11 (v)) and by applying the preceding result of the arithmetic 
mean. □ 

11. Extending means on Hilbert space operators 

We summarize fundamental results of Kudo and Ando [11] (see also [3], [17]) con- 
cerning operator means and their relationship to means on the positive reals. We 
consider continuous means on the positive reals, /i : 1R + x IR + — > IR + , satisfying 

(i) If x < x' and y < y', then fi(x,y) < fi(x',y') (monotonicity) ; 

(ii) fj,(tx,ty) = tfi(x,y) for t,x, y > (homogeneity). 

A homogeneous two-variable function /i can be reduced to a one-variable function 
f(x) = fj,(l,x). This reduction defines a one-to-one correspondence between the 
continuous means satisfying (i) and (ii) and the continuous functions / : R + — > M + 
satisfying 

(0) /(i) = i; 

(1) / is nondecreasing. 

We consider continuous operator means \x on Q, the set of positive operators on a 
Hilbert space, satisfying 

(a) UA<A',B< B', then fi(A,B) < n(A',B') (monotonicity); 

(b) Cn(A,B)C* = n(CAC*,CBC*) for all C invertible (the transformer equal- 
ity), 

where, as usual, A < B means that B — A is positive semidefinite. 

The key result of the theory is that the continuous operator means on Q satisfying 
(a) and (b) are in one-to-one correspondence to the operator monotone functions 
/ : 1R + — > 1R + satisfying (O)-(I), where the correspondence [i f is given by 

H(A,B) = A^fiA'^BA-^A 1 ' 2 . 

Recall that / : 1R + — > 1R + is operator monotone if its extension to Q via the functional 
calculus is monotone. It is this extension that appears in the displayed equality. The 
scalar function / is called the representing function of \i. 

Basic examples of the preceding theory include (i) the arithmetic mean on M + 
with representing function f(x) = (l/2)(l + x) and operator mean (1/2)(A + B), (ii) 
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the geometric mean y/ab on M + with representing function f(x) = y/x and operator 
geometric mean Aj^B = A 1 / 2 (A~ 1 / 2 BA~ 1 / 2 ) 1 / 2 A 1 / 2 , and (iii) the harmonic mean 
2/(a~ 1 + b^ 1 ) with representing function f(x) = 2x/(l + x) and operator harmonic 
mean 2(A~ 1 + S" 1 )" 1 . 

We recall from Section 6 and from [TT] the notion of the iterated composition a * r 
of two operator means a and r on Q. Starting with o~i — a and T\ = r, we define 
inductively the sequences of means {cr n } and {r n } by 

a n+1 {A, B) = a(a n (A, B), r n (A, B)), r n+1 (A, B) = r(a n (A, b), r n {A, B)). 

By Theorem 6.2 of [11] if a and r are continuous, monotonic means satisfying the 
transformer equality and if at least one is neither the left- nor the right-trivial mean, 
then {cr n } and {r n } converge to the same limiting mean (Too, which means that for all 
A, B e f2, lim n a n (A, B) = a^A, B) = lim n r n (A, B), where the limit is taken in the 
weak operator topology. 

Theorem 11.1. Let Q denote the set of positive operators on a Hilbert space, let 
X(A,B) = Aj^B denote the geometric mean of A,B, and let v be continuous, mono- 
tonic mean on Q satisfying the transformer equality that is also coordinatewise p n - 
contractive for < p n < 1 and nonexpansive on the order interval [(l/n)I, nl] for each 
n with respect to the Thompson metric. Then the iterated composition fi = X*u (resp. 
skewed iterated composition) exists, is a coordinatewise p n - contractive, nonexpansive 
mean when restricted to [(l/n)I, nl] for each n, and hence inductively ^-converges to 
a (3 -extension for each n > 2. 

Proof. We consider some fixed order interval T n = [(l/n)I, nl] and [(l/n)J, nl]; 

note that [(l/n)l,nl] is closed and convex with respect to any monotonic mean, in 
particular with respect to A and v. By Proposition 16.41 there exists an iterated com- 
position p n = A|r n * v\v n that is nonexpansive and coordinatewise p n -contractive. 
Clearly if m < n, then p n is an extension of p m . Thus there exists a unique mean p 
that extends all of them. Since any A, B belongs to the domain of some p n , p is the 
iterated composition A * v. 

Since p n is nonexpansive and coordinatewise p n -contractive on T n , it inductively 
admits a /3-extension for each n > 2. It then follows from Corollary 17.81 that a 
/3-extension of p exists inductively for each n > 2. 
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The case of the skewed iterated composition is analogous. □ 

The iterated composition of the operator arithmetic mean and operator geomet- 
ric mean yields the arithmetic-geometric operator mean [8]. Similarly we have the 
harmonic-geometric operator mean. 

Definition 11.2. For two positive definite operators A, B e Q on the Hilbert space E, 
we define the arithmetic- geometric mean or Gauss mean AGM(A, B) to be the iter- 
ated composition of the geometric and arithmetic means, that is, the limit lim n X n (A, B) = 
lim n v n (A, B) , where we define Xi(A, B) = A#B, the geometric mean, U\(A,B) = 
(A + B)/2, the arithmetic mean, and inductively X n+ i(A, B) = X n (A, B)#v n (A, B) 
and u n+1 (A,B) = (l/2)(A n (A, B) + u n {A, B)). 

It is standard that the preceding iteration defining the arithmetic-geometric mean 
of two positive operators converges in the weak operator topology and agrees with 
the one arising from the representing function of arithmetic-geometric mean on the 
positive real numbers [TTJ Section 6]. 

Corollary 11.3. On any order interval [(l/n)I ,nl] the arithmetic- geometric mean 
(resp. the harmonic- geometric mean) is nonexpansive and coordinatewise p n - contractive 
for some p n , < p n < 1. Hence on Q the arithmetic- geometric mean (resp. the 
harmonic-geometric mean) inductively ^-extends to an n-mean for each n > 2. 

Proof. Fix some positive integer m. By Lemma 110.31 the arithmetic mean is nonex- 
pansive and coordinatewise p m -contractive on [(l/m)l,ml] for some < p m < 1. We 
have already remarked that the geometric mean is a convex mean with respect to the 
Thompson metric. The result now follows from Theorem 111.11 □ 

We briefly recall the operator logarithmic mean as discussed in [17]. The loga- 
rithmic mean is defined on IR + by L(a,b) = (b — a)/(logb — log a). Its representing 
function is f(x) = L(l,x) = (x — l)/logx, which is an operator monotone function. 
Hence there exists a corresponding operator logarithmic mean. B. C. Carlson [6] has 
shown that the logarithmic mean on IR + is the skewed iterated composition of the 
geometric and arithmetic mean. It then follows from the theory of operator means 
as developed in Section 6 of |llj . particularly Lemma 6.1 and Theorem 6.2, that the 
operator logarithmic mean is the corresponding skewed iterated composition of the 
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operator geometric and arithmetic means on Q, the set of positive operators on a 
Hilbert space, where the limits are taken in the weak operator topology. However, 
the same arguments applied in the previous corollary to the arithmetic-geometric 
mean viewed as the iterated composition of the geometric and arithmetic means ap- 
ply equally well to logarithmic mean viewed as the skewed iterated composition of 
the geometric and arithmetic mean. We thus obtain analogously the following 

Corollary 11.4. On any order interval [{l/n)l,nl\ the logarithmic mean is nonex- 
pansive and coordinatewise p n - contractive for some p n , < p n < 1. Hence on Q the 
logarithmic mean inductively (3-extends to an n-mean for each n > 2. 

The preceding corollary provides a positive solution to a problem raised by Petz 
and Temesi (p2], [16]) as to whether the logarithmic mean /3-converges and hence 
admits higher dimensional extensions. 

It is easy to obtain the order relation L(A, B) < AGM (A, B) between the logarith- 
mic mean and the arithmetic-geometric mean, which are monotone (Definition 9.1). 
Applying Theorem 19.21 we have the following 

Corollary 11.5. The order relation L n (Ai, . . . , A n ) < AGM n (Ai, . . . , A n ) holds for 
the extended logarithmic and arithmetic- geometric n-means for each n > 2. 

Remark 11.6. In a similar way one can show that 2-means \x and v that satisfy the 
inequality p < v and that can both be derived by some iteration or skew iteration 
of the arithmetic and geometric resp. harmonic and geometric means satisfy p n < v n 
for all n > 2. In this way, for example, one derives the principal results of [10J as a 
corollary to our preceding results. 

We remark in closing that a number of ideas in this paper can be carried over to 
the study of means on the set of positive elements of a C*-algebra, particularly by 
viewing the C*-algebra as a closed subalgebra of the algebra of bounded operators 
on a Hilbert space. For example, one could define the logarithmic mean to be the 
skewed iterated composition of the geometric and arithmetic means and show that it 
inductively /3-extends to all higher dimensions. 



A GENERAL FRAMEWORK FOR EXTENDING MEANS TO HIGHER ORDERS 



35 



References 

[I] T. Ando, C.-K. Li, and R. Mathias, Geometric means, Linear Algebra and Appl. 385(2004), 
305-334. 

[2] W. Ballmann, Lectures on Spaces of Nonpositive Curvature, DMV Seminar Band 25, Birkhauser, 
Berlin, 1995. 

[3] R. Bhatia, Matrix Analysis, Springer- Verlag, New York, 1996. 

[4] J. Borwein and P. Borwein, Pi and the ACM, John Wiley and Sons, 1987. 

[5] M. Bridson and A. Hacfligcr, Metric Spaces of Non-Positive Curvature, Springer, Berlin, 1999. 

[6] B. Carson, The logarithmic mean, Amer. Math. Monthly 79 (1972), 615-618. 

[7] G. Corach, H. Porta, and L. Recht, Convexity of the geodesic distance on spaces of positive 

operators, Illinois J. Math 38 (1994), 87-94. 
[8] J. Fujii, Arithmetico-geometric mean of operators, Math. Japon. 6 (1979), 667-669. 
[9] A. Horwitz, Invariant means, J. Math. Anal Appl. 270 (2002), 499-518. 

[10] Y. Hu, X. Zhang, and Z. Yang, Mixed mean inequalities for several positive definite matrices, 
Linear Algebra Appl. 395 (2005), 247-263. 

[II] F. Kudo and T. Ando, Means of positive linear operators, Math. Ann. 246 (1980), 205-224. 
[12] J. Lawson and Y. Lim, The geometric mean, matrices, metrics, and more, Amer. Math. Monthly 

108 (2001), 797-812. 
[13] — , Symmetric spaces with convex metrics, to appear Forum Math. 
[14] — , Metric convexity of symmetric cones, submitted. 

[15] — , Geometric means and reflection quasigroups, Quasigroups and Related Systems 14 (2006), 
43-59. 

[16] D. Petz, Means of positive matrices: Geometry and a conjecture, Annales Mathcmaticae et 

Informaticae 32(2005), 129-139. 
[17] D. Petz and R. Temesi, Means of positive numbers and matrices, SIAM J. Matrix Anal. Appl. 

27 (2005), 712-720. 

[18] A. C. Thompson, On certain contraction mappings in a partially ordered vector space, Proc. 
Amer. Math. Soc, 14 (1963), 438-443. 

Department of Mathematics, Louisiana State University, Baton Rouge, LA70803, 
USA 

E-mail address: lawsonOmath . lsu . edu 

Department of Mathematics, Kyungpook National University, Taegu 702-701, Ko- 
rea 

E-mail address: ylim@knu . ac . kr 



